Coupling atomistic and continuum hydrodynamics through a mesoscopic model: 

application to liquid water 
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We have conducted a triple-scale simulation of liquid water by concurrently coupling atomistic, 
mesoscopic, and continuum models of the liquid. The presented triple-scale hydrodynamic solver 
for molecular liquids enables the insertion of large molecules into the atomistic domain through 
a mesoscopic region. We show that the triple-scale scheme is robust against the details of the 
mesoscopic model owing to the conservation of linear momentum by the adaptive resolution forces. 
Our multiscale approach is designed for molecular simulations of open domains with relatively large 
molecules, either in the grand canonical ensemble or under non-equilibrium conditions. 



I. INTRODUCTION 

Many relevant properties of condensed matter require 
understanding how the physics at the nanoscale (nm and 
ns) builds up or intertwines with structures and processes 
on the microscale (/im and /zs) and beyond. The so 
called multiscale modeling techniques have been rapidly 
evolving during the last decade to bridge this gap. The 
"multiple-scale" problem is common to many different 
disciplines, and a variety of multiscale models are being 
designed to tackle different scenarios either in solids 1,2 
or soft matter— The main objective of multiscale 

modeling of complex fluids is to study the effect of large 
and slow flow scales on the structure and dynamics of 
complex molecules (e.g. polymers, proteins), or complex 
interactions (e.g. liquid-solid interfaces, wetting fronts, 
structure formation, etc.). In this contexts, multiscale 
modeling is usually based on domain decomposition: a 
small part of the system [O(10nm)] is solved using fully 
fledged (classical mechanics) atomistic detail and it is 
coupled to a (much larger) outer domain, described by a 
coarse-grained (either particle or continuum) model. The 
central idea of these "dual-scale" methods is to solve large 
and slow processes using a computationally low demand- 
ing description, while retaining an atomistic detail only 
where necessary. 

Another important application of domain decomposi- 
tion is the study of open systems. Indeed, in general, 
a dynamic coupling based on domain decomposition re- 
quires to be able to "open up" a molecular dynamics 
(MD) region, in the sense that mass, momentum and 
energy should be exchanged with the exterior in a phys- 
ically sound way. A formulation for flux exchange across 
open boundaries in particle systems is already available 7 , 
and was shown to allow for MD simulations in differ- 
ent type of thermodynamic ensembles. The triple-scale 
scheme presented in this work is equipped with this idea, 
which permits to study the dynamics of confined (yet 
open) molecular systems, which evolve towards the grand 
canonical equilibrium ensemble (see e.gi»). In passing we 
note that the existing Monte Carle 9 (MC) or hybrid MC- 
MD algorithms^ for the grand canonical ensemble can 



only provide restricted dynamical information of the sys- 
tem. 

The first class of schemes based on domain decompo- 
sition to appear were based on particle-continuum cou- 
pling (see^ for a review). In particular, unsteady flow 
can be solved by hybrids based on exchanging the mo- 
mentum flux across the interface (H) between a molecular 
dynamics (MD) domain and a continuum fluid dynamics 
(CFD) solver. One of these schemes (HybridMD) imple- 
ments the open-boundary formulation 7 , which allows to 
open up the MD domain and include mass and energy ex- 
changes across the MD-CFD interface (H). This requires 
consideration of the proper hydrodynamic fluctuations 
at the CFD domai n 11 ! 12 . However, the particle inser- 
tion used by original "open MD" scheme was restricted 
to small solvent molecules, such as argon or water—, due 
to the large steric hindrance of any atomistic complex 
molecule description. 

More recently, another type of domain decomposi- 
tion, based on particle-particle coupling, was developed. 
The Adaptive Resolution Scheme (AdResS)^^ couples 
a coarse-grained particle model with its corresponding 
atomistic description. To do so, the number of degrees of 
freedom of the molecules is adapted (reduced/increased) 
as molecules move across a "transition" layer where the 
all-atom explicit model (ex) and the coarse-grained (eg) 
model are gradually switched on/off, through a hybrid 
model (hyb). The great benefit of AdResS resides in 
making feasible the gradual (on-the-fly) transition of a 
complex molecule description: from a coarse-grained po- 
tential with soft intermolecular interactions to an atom- 
istic one, with the whole set of hard-cores. 

We realized that taking the advantages of HybridMD 
and AdResS should then have a symbiotic effect, poten- 
tially solving most of the limitations of each method. In 
a recent article^ we started to explore in such direction 
and applied the combined AdResS-HybridMD model to 
a liquid of simple tetrahedral molecules. By perform- 
ing molecule insertion within the eg domain, the com- 
bined scheme enables to simulate an open MD system 
and couple its dynamics to a continuum flow descrip- 
tion of the outer region. However, the coupling strategy 
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used in Refi^ 6 . does not avoid some drawbacks already 
present in the original setup of the AdResS scheme^ 5 - 
(methodological improvements to overcome this limita- 
tions are underway). In particular, precise mapping of 
structural and dynamical properties of the eg and hyb 
molecule o 17 i 18 ' 19 i 20 were still required. Hence, any simu- 
lation exploring a new thermodynamic state require new 
calibrations of the coarse-grained models (in practice, 
each simulation is restricted to sample one single ther- 
modynamic state). 

In the present work we show that the coupling geom- 
etry can be modified to yield a more flexible and robust 
AdResS-HybridMD scheme. This new implementation 
avoids the burden associated to the fine tuning of coarse- 
grained layers, thus relieving a great deal of the speci- 
ficity of the coarse-grained model. This is important 
not only from the computational standpoint, but also 
because it opens a route to consider processes following 
a thermodynamic path. A general completion of this 
route requires the inclusion of the energy exchange into 
the combined scheme, and some possible solutions are 
hereby suggested. Section HT1 briefly introduces the (dual- 
scale) hybrid models (HybridMD and AdResS). Coupling 
strategies are explained in Sec. (jllip and simulations and 
results are presented in Sees. (|IV|) and (|v). Some con- 
clusions are given in (|VI[) . 



II. HYBRID MODELS 
A. Particle-continuum hybrid (HybridMD) 

The HybridM D 11 ! 12 scheme couples the hydrodynamic 
of a particle region, here called "molecular dynamics" 
(MD) domain, with a continuum fluid dynamic descrip- 
tion (CFD) of the external fluid. The essential quantity 
exchanged between CFD and MD is the momentum flux 
across the MD-CFD interface H (see Fig. I); which can 
be casted as Jh • n, where 3h is the local pressure tensor 
and n is the unit vector normal to the H-interface, whose 
area is A. The momentum flux is transferred to the MD 
domain by imposing an external force F ext = AJh • n to 
a particle buffer (the overlapping domain B in Fig. I) ad- 
jacent to the MD domain. Molecules are free to cross the 
H-interface, from or towards the buffer, but once in B, 
each molecule i feels an external "hydrodynamic" force 
distributed according to Ff xt — g(xi)F ext / Y^ieB 9{ x i)i 
where x is the coordinate normal to H. Several options 
can be chosen for the distribution function 2 ^ we used a 
step function g(x) = Q(x — x Q ), as in Ref^ 2 -. 

Through the interface H, the CFD domain receives ex- 
actly the same amount of momentum as the MD system 
does, but in opposite direction, thus ensuring conserva- 
tion. Also, the particle's mass crossing H is injected into 
the CFD domain, via a relaxation procedure 11 . This 
means that conservation of total mass and momentum 
only apply to the system MD+CFD. In other words, 
strictly speaking B is not part of the total system^ and 



can be thought of as a particle reservoir where we apply 
flux boundary conditions to the MD system. Finally, we 
note that the only "microscopic" information required at 
the CFD level is the equation of state and viscosities of 
the atomistic fluid (Here, we assume that the constitutive 
relation is known). 

B. Adaptive Resolution Scheme (AdResS) 

As stated, the AdResS scheme couples an atomistic do- 
main and a coarse-grained particle region by adapting, 
on-the-fly, the molecular description of those molecules 
moving across both domains. This is clearly illustrated 
in Fig. fb,c, where the atomistic domain is labeled as 
ex, the coarse-grained region as eg and hyb denotes the 
transition region between both domains. In this transi- 
tion regime, the force acting on a molecule is an hybrid 
of the explicit and coarse-grained forces, 

F aP = w{X a )w(X )Y a $ m + [l-w{X a )w{Xp)}¥ c $, (1) 

where F„^ is the intermolecular force acting between cen- 
ters of mass of molecules a and ft, placed at x = X a 
and Xj3 in the coupling coordinate, F^ m is the sum 
of all pair intermolecular atom interactions between ex- 
plicit atoms of the molecules a and ft, FJj™ = —VapU ™ 
is the corresponding coarse-grained force acting on the 
center of mass of the molecule (cm), and w is the weight- 
ing function determining the degree of resolution of the 
molecules. The value w = 1 corresponds to the ex region, 
w = to the eg region, while values < w < 1 corre- 
sponds to hybrid (hyb) models. In this study we used the 
same functional form of w as in RefJ£. 

Each time a molecule leaves (or enters) the coarse- 
grained region it gradually gains (or loses) its equilibrated 
vibrational and rotational DoFs while retaining its lin- 
ear momentum. Note that the change in resolution car- 
ried out by AdResS is not time-reversible as a given eg 
molecule corresponds to many orientations and configu- 
rations of the corresponding ex molecule. Since time re- 
versibility is essential for energy conservation^ 2 -, AdResS 
does not conserve energy. In particular, the force in Eq. 
|T| is in general not conservative in the hyb region (i.e., in 
general § F Q ^ • dr ^ 0) 23 i 24 . Hence, to supply or remove 
the latent heat associated with the switch of resolution 
we employ a standard DPD thermosta t 25 ' 26 acting at the 
eg and hyb regions. Note that the thermostat forces do 
not enter into the AdResS interpolating scheme, Eq. (fT|). 
Instead they are added to the AdResS^ 5 -. 

III. COUPLING STRATEGIES 

The AdResS-HybridMD scheme can be applied in sev- 
eral contexts. In general terms, the scheme solves the 
hydrodynamic coupling of a MD system with external 
(CFD) flow. This goal requires hydrodynamic consis- 
tency (e.g. momentum conservation). Also, a reduced, 
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but yet important, application of the combined scheme 
consists on the study of the equilibrium molecular dy- 
namics of open systems with relatively large molecules. 
For this sake thermodynamic consistency is required: 
in particular sampling the grand-canonical ensemble re- 
quires proper mass fluctuations across the simulation 
boundaries. Confined systems are a relevant example 
of this sort of application (see e.g. Ref. 8 ), for which the 
role of CFD domain can be simplified to just provide the 
external pressure (and temperature) of the external mass 
reservoir. 

In the same way, the geometry used for the AdResS- 
HybridMD coupling allows for a certain flexibility; the 
key issue being the location of the interface H connecting 
the MD and CFD domains (see Fig. 1). In a previ- 
ous worki^ we proposed to place the H interface within 
the coarse-grained domain (see Fig. la). This setup is 
useful to include hydrodynamics in simulations requir- 
ing a particle-based multiscale description (e.g. using 
AdResS). An example of application could be the study 
of the dynamics of ions nearby a charged surface un- 
der flow conditions; here the ex layer would describe 
the physics near and at the surface, while away from it, 
charged eg particles will eventually transfer the external 
flow from the outer CFD region. 

On the other hand, many studies are focused on the 
atomistic region and do not really require a particle-based 
multiscale description around it. In these cases, a more 
logical setup is to place the interface H within the atom- 
istic domain ex (see Fig. lb,c). In this paper we explore 
the benefits of this second coupling strategy which, as 
stated, might open a route to consistently couple the hy- 
drodynamics and thermodynamics of qualitatively differ- 
ent levels of description, 

A. Coarse-grained buffer 

A detailed description of this setup which locates the 
hybrid interface H into the eg domain (see Fig. la), can 
be found in Ref J^; we now briefly summarize its require- 
ments. In this setup the original AdResS scheme is im- 
plemented inside the the MD domain. Thus, the first 
requirements of this coupling geometry are demanded 
by AdResS^. In particular, to guarantee a similar fluid 
structure along the AdResS layers, one needs to calibrate 
the radial distribution function and equation of state of 
the coarse-grained model eg and hyb so as to fit the atom- 
istic fluid values. Extra pressure corrections require also 
the same type of calibrations for the intermediate hybrid 
model w — 1/2. On the other handiS, hydrodynamic 
consistency demands to take care of the viscosities of the 
eg and hyb models, and fit them to the ex fluid valued 
(as shown in Ref^ this does not guarantees a perfect fit 
of the diffusion coefficients). In summary, for each ther- 
modynamic state considered, one needs to perform the 
following calibration steps: 

1. Calibrate the effective potential U^2. of the eg 



model 17 so as to fit the center-of-mass radial dis- 
tribution function (RDF cm ) and pressure of the ex 
model 3 (solid line in the inset of Fig. \2§ ■ 

2. The interface pressure correction^ 7 - is used to sup- 
press density oscillations in the hyb layer. This 
requires the calibration of the effective potential 
C^u=o 5 f° r the hybrid model with w = 1/2 (dashed 
line in the inset of Fig. [5]). 

3. Measure the viscosity of the eg model. Then cali- 
brate the transverse friction coefficient of the trans- 
verse DPD thermostat^ so as to match n cg and rj ex . 
This requires several viscosity calculations^. 



B. Adaptive resolution buffer 

The setup discussed in this work consists on placing the 
hybrid interface H inside the atomistic domain (see Fig. 
lb,c). In this setup AdResS works inside the buffer. The 
MD region is thus purely atomistic and this fact brings 
about an important benefit: around the hybrid inter- 
face H, all relevant fluid properties are properly defined. 
These include the fluid transport coefficients, energy and 
pressure equation of states and the corresponding mass 
and pressure fluctuations. In this context, MD-CFD cou- 
pling by HybridMD is well defined 11 . On the other side, 
by placing AdResS into the particle buffer it behaves as 
an adaptive resolution buffer, where the atomistic DoF's 
are gradually inserted into the MD (atomistic) region of 
interest. Note that the eg layer, with soft interaction po- 
tentials, is placed near the buffer end (see Fig. lb,c), so 
the insertion of complex molecules into the system is still 
an easy task, using standard insertion routines^. 

The second relevant benefit of this coupling strategy 
is that it permits to properly impose the external pres- 
sure and stress into MD, without having to perform a 
fine tuning of the structural and dynamic properties of 
the eg and hyb models. In other words, the AdResS- 
HybridMD scheme does not rely on the specificity of the 
coarse-grained description anymore. The reason for this 
is simple: the AdResS scheme conserves linear momen- 
tum. This means that any external momentum flux will 
be properly transferred across the AdResS layers to the 
atomistic core. As the external "hydrodynamic" force 
used in the HybridMD scheme decouples from any inter- 
molecular interaction (the total force on a molecule at 
the buffer B is, F a = F^f* + J2a^ap), the transfer of 
mechanical variables (pressure and stress) should be ro- 
bust against the details of the eg and hyb models. We 
now present simulations to prove this claim. We will also 
consider any effect on the liquid structure near the in- 
terface H and check for proper mass fluctuations across 
H. 
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IV. SIMULATIONS 

The present hybrid simulations were implemented for 
the flexible TIP3P water model, partly because of the 
relevance of water and also because it has well known 
structural and dynamical propertie o 29 i 30 ' 31 i 32 . In the 
remainder of the paper we use reduced Lennard-Jones 
units corresponding to the Oxygen atom: mass m = 
?7io = 16 a.u., oxygen-oxygen interaction energy e = 
( oo = 0.152073kcal/mol, and diameter a — ao — 
3.1507A Simulations were done at ambient tempera- 
ture, T = 300K, which corresponds to ksT/e = 3.92 
in reduced units. TIP3P-water density is 1.02 gr/cm 3 , 
which corresponds to p = 1.20 m/a 3 . Most particle sim- 
ulations were done in boxes of total (MD+B) volume 
24.5 x 6. 18 x 11. 12c 3 , although in order check for finite size 
effects, boxes with larger dimensions were also used. The 
volume of the MD domain was V = 10.5 x 6.18 x 11.12cr 3 . 
and it contained about 865 water molecules. Long-range 
electrostatic forces are computed using the reaction field 
method, in which all molecules outside a spherical cav- 
ity of a molecular based cutoff radius R c — 2.86cr are 
treated as a dielectric continuum with a dielectric con- 
stant erf — 8 27 i 33 ' 34 . The finite volume solver for the 
CFD domain was feeded with the equation of state for 
flexible TIP3P water reported in Ref? 3 -^. Finally, the mi- 
croscopic part of the stress tensor at the hybrid interface 
H was measured according to the mesoscopic approach 
explained in Ref. 12 ; i.e. by evaluating velocity gradients 
from the cell-averaged particle velocities and using the 
Newtonian constitutive relation. 



V. RESULTS 

Coarse-grained buffers. As expected, results proved 
the correct behavior of the combined scheme, both in 
terms of the structure and hydrodynamics at the MD 
domain. However calibration steps, using the recipe ex- 
plained in Sec. IIII Al involved significant work. Strong 
corrections of the eg and hyb viscosities^ were required 
(without transverse DPD thermostat, r\ cg is about five 
times smaller than r\ ex ). In these calibration steps we 
used HybridMD as a rheometer— . 

Adaptive resolution buffers. In what follows we focus 
on the result obtained for the second, much lighter setup. 
It is important to stress that these simulations were done 
using an effective potential for the eg model which was 
deliberately not accurately fitted to reproduce the all- 
atom RDF cm (see Fig. [2]). Moreover, steps 2 and 3 of 
the protocol of Sec. IIII Al were also avoided, meaning 
that the shear viscosities at the eg and hyb layers result 
in much smaller values than the atomistic (ex) one. 



A. Liquid structure 

We first compare, in Fig. [2 the local structure of the 
liquid inside the MD region of the triple-scale model with 
that obtained from all-atom simulations within periodic 
boundaries. The agreement is perfect, indicating that 
the unfitted liquid structure in the buffer (dashed line in 
Fig. [5]) does not affect the proper liquid structure inside 
the interest (MD) region. 

B. Hydrodynamics 

The hydrodynamic behavior of the triple-scale scheme 
was tested considering both steady and unsteady flows. 
Figure [3] shows the density and velocity profiles in the 
particle region (MD+B) obtained at the steady state of 
a simple Couette flow. The density (and temperature, 
not shown) profile inside the MD region is flat and con- 
firm that the triple-scale scheme furnishes an homoge- 
neous equilibrated liquid bulk. The expected linear ve- 
locity profile inside the MD domain indicates that the 
transverse momentum is correctly transferred across the 
triple-scale fluid model. We also conducted simulations 
of Stokes flow (an oscillatory shear flow driven by the 
oscillatory motion of a wall along its plane direction). 
Figure [4] shows the velocity in one of the MD-cells of the 
system corresponding to a flow with a period 300r. The 
hybrid result perfectly agrees with with the deterministic 
solution of the Navier-Stokes equation (red solid line in 
Fig. [4|. The large viscosity of liquid water induces a fast 
transfer of momenta across the buffers. Hence, even for 
faster shear rates, no trace of phase delay between the 
MD and CFD velocities is observed^, (see the inset of 
Fig. ij. 

C. Mass fluctuations 

One of the interesting properties of the AdResS- 
HybridMD approach is that the MD region becomes an 
open system, which exchanges mass with its surround- 
ings. As stated before, to that end it is quite impor- 
tant to check that mass fluctuation across the MD bor- 
der (H) is consistent with the grand canonical prescrip- 
tion. We measured the mass variance inside the MD 
domain and compared it with the grand canonical (GC) 
result, Var[p] = pksT/ (Vc^), where V is the system's 
volume and c\ = {dP/dp)T is the squared isothermal 
sound velocity (related to the isothermal compressibility, 
(3t = ( c tP) -1 )- Taking the sound velocity for the flexible 
TIP3P water reported in Ref, 31 , c T = 7.38(e/m) 1 / 2 and 
the mass density p — 1.20m/ a 3 (recall that m — mo), 
the GC prediction for V = 3.5 x 6.18 x 11.12cr 3 is Var[p] = 
0.0187. Inside the MD domain we obtained Var[p] = 
0.020 ± 0.002 within different slices of the same volume. 
This is a quite good agreement, considering the smallness 
of V. In a larger volume V — 10.5 x 6.18 x 11.12(T 3 , the 
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triple-scale result Var[p] = 0.011 ± 0.005 is even closer to 
the GC prediction Var[p] = 0.0108. 

D. Energy 

Although in this work we do not solve the energy ex- 
change across H, we are in position to provide some ar- 
guments indicating that this task is solvable using the 
adaptive resolution buffer setup. Energy exchange re- 
quires three properties: 1) energy should be properly de- 
fined across MD, 2) the scheme should allow to change 
the thermodynamic state of the system and 3) one should 
be able to control the amount of energy per unit time in- 
serted into MD. Concerning 1), it is important to stress 
that by placing the AdResS scheme inside the buffer one 
ensures that the energy is a well defined quantity over 
the total (MD+CFD) system. Also 2) is satisfied be- 
cause the adaptive resolution buffer setup does not rely 
on the specificity coarse-grained model: this means that 
it should be possible to change the thermodynamic state 
of the MD domain (e.g. the mean temperature) without 
re-calibrating the eg layer (we have tested this in simu- 
lations using tetrahedral molecules). Finally, we believe 
that 3) is solvable, but it will require a further devel- 
opment of the algorithm. One possible way could be 
to implement the flux-boundary conditions developed in 
Refi We expect to present an algorithm of this sort in 
future works. 



VI. CONCLUSIONS 



To conclude, we have presented a flexible and ro- 
bust hybrid scheme for hydrodynamics of molecular liq- 
uid, which combines atomistic, mesoscopic and contin- 
uum models. This triple-scale scheme uses a flux based 
particle-continuum hybrid to couple the atomistic core 
and continuum sides of the system, while it generalizes 
the role of the particle buffer to allow for a gradual change 
in molecular resolution: from an all-atom at the core to a 
mesoscopic one near the buffer end. Structure and hydro- 
dynamic of the core (MD) region were shown to be robust 
against changes in the choice of the mesoscopic model, 
greatly reducing calibration burdens. Further extensions 
to allow for energy exchange and multiple species will be 
explored in future works. 
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FIG. 1: Coupling strategies for the AdResS-HybridMD 
scheme: (a) Coarse-grained buffer and (b), (c) Adaptive res- 
olution buffer. In each figure, the bottom panel depicts the 
decomposition of the whole system (MD+CFD); where MD 
stands for the molecular dynamics region, surrounded by con- 
tinuum fluid dynamics (CFD) domains solved via the finite 
volume method^. The MD-CFD coupling is solved by the 
HybridMD scheme—, based on the exchange of momentum 
flux across the interface H. Pressure and stress are imposed 
into MD via external forces acting on particles at the buffers 
B. The AdResS scheme^ (see the middle figure) gradually 
adapts the atomic resolution of the molecules: from all-atom 
(ex) to coarse-grained (eg) descriptions, passing through a hy- 
brid (hyb) model. AdResS and HybridMD can be combined 
in two ways depending on location of the interface H: either 
using a Coarse-grained buffer (a), see Ref.— ; or an Adap- 
tive resolution buffer (b) and (c), explored in this work. Fig- 
ure (c) is an illustration of this triple-scale scheme for liquid 
water. The hydrodynamic coupling is made along x direc- 
tion, (finite volume cells are Ax = 3.5a wide) and the system 
is periodic in the orthogonal directions. The atomic resolu- 
tion of water molecules is gradually switched on as they move 
across the buffer, which is la long. The hyb region is 3.5<r 
and its distance to H is about la. A standard DPD thermo- 
stat at T = 300K is used for the eg and hyb domains, with 
a friction constant £ = 0.5m/r. Information between MD 
and CFD is exchanged after every fixed time interval At c , 
with At c = ncFDAtcFD = nMD&t. Here we typically used 
At c = At ~ 0.03r and St = 0.0003r = 0.5fs (small enough to 
recover O-H vibrational motion). 
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FIG. 2: Comparison between RDF cm s obtained in a all- 
atom (ex model) periodic boundary simulation of the flexible 
TIP3P water at ambient conditions and within the MD region 
of a triple-scale simulation. At the buffer we used a eg model 
whose RDF cm (dashed line) was not fitted to the all-atom re- 
sult. The inset shows the effective potentials and t/^,=o.5 
used in the first protocol (see text), which correctly reproduce 
the all-atom RDF cm : two minima corresponding to the first 
and second hydration shells of the liquid are observed. 




20 25 
position (a) 



FIG. 3: Density profile (solid line) and velocity distribution 
(squares) across the particle domain in a steady Couette flow 
(the dashed line is the expected linear profile across the MD 
region) . 
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FIG. 4: Time evolution of the velocity at one finite volume 
cell within the MD region in an oscillatory shear flow. For 
comparison, the deterministic Navier- Stokes solution is shown 
in red lines. The inset shows velocity in one MD cell at the 
start-up of a Couette flow. 
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